Mpemba effect and phase transitions in the adiabatic cooling of water before freezing 
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An accurate experimental investigation on the Mpemba effect (that is, the freezing of initially 
hot water before cold one) is carried out, showing that in the adiabatic cooling of water a relevant 
role is played by supercooling as well as by phase transitions taking place at 6 ± 1°C, 3.5 ± 0.5°C 
and 1.3 ± 0.6°C, respectively. The last transition, occurring with a non negligible probability of 
0.21, has not been detected earlier. Supported by the experimental results achieved, a thorough 
theoretical analysis of supercooling and such phase transitions, which are interpreted in terms of 
different ordering of clusters of molecules in water, is given. 



A well-known phenomenon such as that of the freezing of 
water has attracted much interest in recent times due to 
some counter-intuitive experimental results [l[ and the 
apparent lacking of a generally accepted physical inter- 
pretation of them 0, 0], [H ■ These results consist 
in the fact that, many times, initially hot water freezes 
more quickly than initially cold one, a phenomenon which 
is now referred to as the Mpemba effect (for a short his- 
torical and scientific survey see the references in |3j). The 
observations sound counter-intuitive when adopting the 
naive, simple view according to which initially hot water 
has first to cool down to the temperature of the initially 
cold one, and then closely follow the cooling curve of 
the last one. The effect takes place even for not pure 
water, with solutions or different liquids (the original 
Mpemba observation occurred when he tried to make an 
ice cream). 

Several possible physical phenomena, aimed to explain 
such observations, have been proposed in the literature, 
mainly pointing out that some change in water should 
occur when heated Q H- 

However, such explanations cannot be applied if some 
precautions are taken during the experiments (whilst the 
Mpemba effect has been observed even in these cases) 
and, in any case, calculations do not seem to support 
quantitatively the appearance of the effect (see references 
in |). 

Some novel light has been introduced in the discussion, 
in our opinion, in Ref. 0], where the Mpemba effect 
has been related to the occurrence of supercooling both 
in preheated and in non-preheated water. Initially hot 
water seems to supercool to a higher local temperature 
than cold water, thus spontaneously freezing earlier. As 
a consequence, such a scenario, apparently supported by 
experimental investigations, points toward a statistical 
explanation of the effect, neither the time elapsed nor 
the effective freezing temperature being predictable. 

Here, we prefer to face the problem by starting from 
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what is known about the freezing process, rather than 
the cooling one. 

In general it is known that, for given values of the ther- 
modynamic quantities (for example the volume and the 
energy), a physical system may exist in a state in which 
it is not homogeneous, but it breaks into two or more 
homogeneous parts in mutual equilibrium between them. 
This happens when stability conditions are not fulfilled, 
so that a phase transition occurs; it is, for example, just 
the case of water that, at the pressure p of 1 atm and at 
temperature T of 0°C, becomes unstable. 

When liquid water is cooled, the average velocities of 
its molecules decreases but, even if the temperature goes 
down to 0°C (the fixed temperature where liquid and 
solid phases coexist) or lower, this is not a sufficient con- 
dition for freezing to start. In fact, in order that ice 
begins to form, first of all some molecules of the liquid 
water should arrange in a well-defined order to form a 
minimum crystal and this, in the liquid state, may hap- 
pen only randomly. Second, such starting nucleus has 
to attract further molecules in the characteristic loca- 
tions of the crystalline structure of ice, by means of the 
interaction forces of the nucleus with the non-ordered 
molecules in the liquid. Nucleation and crystal growth 
processes are both favored at temperatures lower than 
0°C, so that supercooling of liquid water is generally re- 
quired before its effective freezing. In fact, in pure water, 
only molecules in the liquid with statistically lower veloci- 
ties can arrange the initial nucleus and, furthermore, only 
slow moving molecules are able to join that cluster and 
put their kinetic energies into potential energy of bond 
formation. When ice begins to form, these molecules are 
removed from those attaining to the given Maxwell dis- 
tribution for the liquid water, so that the average speed 
becomes larger, and the temperature of the system rises 
to 0°C (obviously, the temperature is set at the value 
where the continuing exchange of molecules is equal in 
terms of those joining and those leaving the formed crys- 
tal surface). 

Thus supercooling is, de facto, a key ingredient in the 
freezing process, although supercooled water exists in a 
state of precarious equilibrium (water is in a metastable 
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state). Minor perturbations such as impurities or other 
can trigger the sudden appearance of the stable crys- 
talline phase for the whole liquid mass, again with the 
release of the entire crystallization heat (melting heat) 
which increases the temperature of the freezing liquid to 
the normal 0°C one. 

In general, when a system is in a metastable state, 
sooner or later it will pass to another stable state. In 
water, density and entropy fluctuations favor the forma- 
tion of crystallization nuclei but, if the liquid constitutes 
a stable state, such nuclei are always unstable and will 
disappear with time being. However, since the fluctua- 
tions become more pronounced the lower the tempera- 
ture, if water is supercooled, for sufficiently large nuclei 
they will result to be stable and grew with time, becom- 
ing freezing centers. The starting of the phase transition 
is thus determined by the probability of appearance of 
those nuclei, and the reported Mpemba effect could be 
simply a manifestation of this process. 

We have calculated just this probability P as function 
of the absolute temperature T of the metastable phase 
(the one at which the nucleus is in equilibrium with the 
liquid), obtaining the following result llj : 



P = — exp 



{T-T, 



(1) 



Here is the equilibrium temperature of the liquid-solid 
phase, a is a dimensionless normalization factor and /3 is 
a constant whose expression is given by 



16ttt 3 v 2 



(2) 



where r is the surface tension, v the molecular volume 
of the crystallization nucleus, Q the molecular heat of 
the transition from the metastable phase to the nucleus 
phase, and k is the Boltzmann constant. Just to give 
an idea of the macroscopic value of the constant (3, let 
us note that t 3 v 2 = W s 3 urf is the cube of the work done 
by the surface forces, and by assuming that Q ~ kT* we 
may write: 



167T /Wsurf 

3 ^ Q 



(3) 



that is the constant (3 is ruled by the ratio W SUT {/Q. 

The probability P has a minimum at the liquid-solid 
equilibrium temperature T» and increases for decreasing 
temperature, as expected. From the formulae above it 
is clear that the probability for nucleation, and thus the 
onset of the freezing process as well, is enhanced if the 
work done by the surface forces (or the surface tension 
itself) is lowered in some way. In normal daily conditions 
when a commercial refrigerator is employed, this is eas- 
ily induced in two simple ways: either by the presence 
of impurities, when solutions (such as an ice cream so- 
lution, as in the Mpemba case) are used as the freezing 





V= 20 cm 3 


V= 50 cm 3 


V= 65 cm 3 


V= 80 cm 3 


Psc 


1 


0.28 





0.46 




T c = -8°C 


T c = -14°C 


T c = -22°C 


T c = -26°C 


Psc 


0.75 


0.50 





0.11 



TABLE I: Probabilities for the occurrence of supercooling for 
different volumes V of the sample and different temperatures 
of the cryostat T c . 



liquid instead of pure water, or by fluctuations of the ex- 
ternal pressure or temperature, caused in the commercial 
refrigerator itself. This explains why no appreciable su- 
percooling is observed in normal situations. Obviously, 
the most direct way to induce freezing in supercooled 
water is to introduce an external body in it, in order to 
directly lower the surface tension. 

We have thus performed an accurate experimental in- 
vestigation, accounting for a total of about one hundred 
runs, aimed to clarify the phenomenology of the Mpemba 
effect and its interpretation. In the first part of our ex- 
periments we have tested all the above qualitative pre- 
dictions about supercooling, by studying the cooling and 
freezing of tens of cm 3 of normal water in a commer- 
cial refrigerator, in daily operation conditions. The key 
point, in fact, is not to obtain the most favorable physical 
conditions, employing sophisticated setups, but rather 
to reproduce the Mpemba conditions, that is adiabatic 
cooling (with commercial refrigerators) of not extremely 
small quantities of water. We have used an Onofri re- 
frigerator for the cooling of double distilled water and a 
NiCr-Ni thermocouple as a temperature sensor (Leybold 
666193), interfaced with a Cassy Lab software for data 
acquisition. 

For fixed temperatures of the cryostat we have indeed 
observed supercooling in our samples, with the freezing 
occurring just along the lines predicted above. In par- 
ticular, during the supercooling phase we have induced 
a number of small perturbations in our samples, namely, 
variations of external pressure or temperature, mechan- 
ical perturbations or introduction of an external macro- 
scopic body (a glass thermometer held at the same tem- 
perature of the sample). In all these cases we have regis- 
tered the sudden interruption of the supercooling phase 
and a practically instantaneous increase of the temper- 
ature to the value of 0°C, denoting the starting of the 
freezing process. Conversely, if no perturbation is in- 
duced (or takes places) the water reached an equilibrium 
with the cryostat at temperatures up to about — 30°C 
(lasting also for several thousands of seconds). 

We have then verified that when the freezing process 
started from the supercooling phase, the Mpemba effect 
took place with a probability in agreement with that re- 
ported in Ref. 0]. 

In about half (with a total probability of 0.47) of the 
runs performed we have detected a supercooling phase. 
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FIG. 1: Left: Cooling curves for V = 20 cm 3 and T c = -8°C. 
Right: The fitted time duration of the phase transition at 
3.5°C as function of the volume V of the samples, for different 
temperatures T c of the cryostat. 

In Table U we report the observed probability Psc for 
the occurrence of supercooling for different volumes V of 
the water sample and for different temperatures T c of the 
cryostat. We find the data to be fitted by a straight line, 
denoting (in the range considered) a linear decreasing 
of Psc for decreasing temperatures of the cryostat and 
for increasing volumes of the samples, this probability 
reaching the maximum Psc — 1 for T c — 0°C (and V = 

o). 

An interesting feature of what we have observed is the 
sensible appearance of iced water in our samples. In fact, 
when supercooling did not occur, the ice started to form 
around the walls of the beaker, while the inner parts were 
still in a liquid form, as usually expected. Instead the im- 
mediate freezing of supercooled water involved the whole 
sample, this showing a very peculiar symmetric form. 
We have used cylindrical beakers with the temperature 
sensor in their periphery, near the walls; the observed 
structure was a pure radial (planar) one, with no liquid 
water and radial filaments of ice from the center of the 
beakers to the walls (in one case we have been also able 
to take a low resolution picture of this, before its destruc- 
tion outside the refrigerator). 

However, although supercooling plays a relevant role in 
the manifestation of the Mpemba effect, the things are 
made more complicated by the occurrence of other statis- 
tical effects before the temperature of the water reaches 
the value of 0°C. This comes out when an accurate mea- 
surement of the cooling curves is performed (some exam- 
ples of what we have obtained during the second part of 
our experiments are reported in Fig. [1]). 

According to a simple naive model, the heat exchange 
from the water sample (at initial temperature To) to the 
cryostat (at fixed temperature T c ) is described by the 
equation 

CdT = 6(T C -T )dt, (4) 

where C and 5 are the thermal capacity and the heat 
conductivity of the water, respectively. Thus by solving 
the differential equation in Q, the following expression 



T c = -8 ± 2°C 





V= 20 cm 3 


V= 50 cm 3 


V= 65 cm 3 


V= 80 cm 3 


Ati (s) 
At 2 (s) 
At 3 (s) 


7±1 
11±6 
12 ±6 


220 ± 100 


500 ± 170 
70 ±30 


630 ± 160 


T c = -14±2°C 




V= 20 cm 3 


V= 50 cm 3 


V= 65 cm 3 


V= 80 cm 3 


Ati (s) 
At 2 (s) 
At 3 (s) 


8±3 
7±4 


130 ± 80 


480 ± 160 


37 ±1 
500 ± 60 


T c = -22 ± 1°C 




V= 20 cm 3 


V= 50 cm 3 


V= 65 cm 3 


V= 80 cm 3 


Ati (s) 
At 2 (s) 
At 3 (s) 




63 ±1 
170 ± 100 




7± 1 
130 ± 70 


T c = -26 ± 1°C 




V= 20 cm 3 


V= 50 cm 3 


V= 65 cm 3 


V= 80 cm 3 


Ati (s) 
At 2 (s) 
At 3 (s) 




3.5 ±0.7 

3± 1 
200 ± 70 


320 ± 70 
1±1 


210 ± 170 



TABLE II: Time duration of the phase transitions at 6°C 
(Ati), 3.5°C (Ati) and 1.3°C (Ats) for different volumes V 
of the sample and different temperatures of the cryostat T c . 

for the temperature as function of time t is obtained: 

T = T c -(T c -T )e-*/ T , (5) 

where r = C/S is a time constant measuring the cool- 
ing rate of the sample. However, although the overall 
dependence of T on time is that expressed by Eq. |5|), 
our experimental data clearly reveal the presence of three 
transition points before freezing (or supercooling), where 
r changes its value. This transitions occur at tempera- 
tures Ti = 6±1°C, T 2 = 3.5±0.5°C and T 3 = 1.3±0.6°C 
with a probability of Pi = 0.11, P 2 = 0.84 and P 3 = 0.21, 
respectively. The time duration At of each phase tran- 
sition, during which the temperature keeps practically 
constant [T^, depends on the volume of the sample and 
on the temperature of the cryostat. The data we have 
collected are summarized in Table |TTJ For the phase 
transition at T 2 these data show a linear increase of 
At 2 with T c and a quadratic one with V; in Fig. [1] we 
give the fitting curves corresponding to best fit function 
At 2 = (a+bT c )V 2 . Instead, for the other two phase tran- 
sitions no sufficient data are available in order to draw 
any definite conclusion on the dependence on V and T c , 
though Ati and At 3 appear to be shorter than At 2 . 

The occurrence of these phase transitions is likely re- 
lated to the formation of more or less ordered structures 
in water, resulting from the competition between long- 
range density ordering and local bond ordering maxi- 
mizing the number of local bonds Q. The anomalous 
density maximum at about 4°C (which we observe here 
at T 2 — 3.5 ± 0.5°C) is, for example, explained just in 
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T = — 8°C 




T c = — 22°C 


T c = — 26°C 


n 


(8) 


600 ± 110 


680 ± 100 


1000 ± 110 


950 ± 190 


T-2 


00 


1080 ± 260 


1060 ±170 


530 ± 90 


570 ±3 


r-.i 


00 


1590 ± 930 


1520 ± 730 


270 ± 50 


220 ± 80 


T4 


(b) 


620 ± 480 


500 ± 180 


150 ± 30 


640 ± 490 



TABLE III: Time constants n (T < Ti), r 2 (Ti < T < T 2 ), 
t 3 (T 2 < T < T 3 ), r 4 (T > T 3 ) of the cooling curves before 
and after the three phase transitions detected, for different 
temperatures of the cryostat T c . 



term of this: as water is cooled, the local specific volume 
increases due to the progressive increase in tetrahcdral 
order, so that the entropy, that always decreases upon 
cooling, at 4°C becomes anticorrelated with the volume, 
resulting in an inversion (from positive to negative) of 
the thermal expansion coefficient and a corresponding 
density maximum Similar explanations in terms of 
different ordering could apply also to the other two tran- 
sitions we have observed, but an exhaustive discussion 
of them, which would require more experimental data, 
is beyond the scope of this Letter. We only note that, 
while the first transition at Ti = 6 ± 1°C seems related 
to the effect observed in Ref. at 8°C, to the best of 
our knowledge no other author has reported the one at 
T3 = 1.3 ± 0.6°C (which, as mentioned, occurs with an 
appreciable probability of 0.21). 

The observed mean values of the four time constants 
of the cooling curves, before and after the three phase 
transitions, are reported in Table IIIII for different val- 
ues of T c . All the time constants are approximately in- 
dependent on the volume V, in disagreement with the 
naive model discussed above which predicts an increase 
of t with the thermal capacity. Instead they depend lin- 
early on T c , showing a negative slope for t\ and positive 
ones for T2,T3,T4 and a finite value for T c = 0°C. Note 
that (in the naive model) the ratios of the different time 
constants, at fixed volumes, give the (inverse) ratios of 
the heat conductivities in the different ordered phases 
(all these ratios decrease with the cryostat temperature), 
which are directly related to microscopic quantities like 
the size and average velocity of the ordered clusters of 
molecules in water. 

Coming back to the Mpemba effect, it is easy to see 
that Eq. {H} predicts that, for constant t, initially hot 
water reaches the freezing point later than initially cold 
water. However, from what just discussed, in general this 
could be no longer true if the time constant changes its 
value during the cooling process (the slope of the cooling 
curves changes), or phase transitions before freezing oc- 
cur (with time durations sufficiently long/short). In ad- 
dition to these effects, the reaching of the freezing point 
does not automatically guarantees the effective starting 
of the freezing process, since relevant supercooling may 
take place, thus statistically causing the freezing of ini- 



tially hot water before cold one. 

From the data we have collected we have verified that, 
for given V and T c , in many cases no inversion between 
the cooling curves happens before the freezing point, ir- 
respective of the change in the value of r or the time 
duration of the phase transitions. Nevertheless we have 
as well realized that this is mainly due to the not very 
large difference between the initial temperatures of the 
samples, and in few cases (among those studied by our- 
selves) it cannot be applied, the largest effect causing the 
inversion being the phase transition at Ti- 

In conclusion our experimental results, and their in- 
terpretation reported here, clearly point out the statisti- 
cal nature of the Mpemba effect (as already realized in 
Q), whose explanation is given in terms of transitions 
between differently ordered phases in water and super- 
cooling. The very detection of such phenomena seems 
to require the cooling to be adiabatic (as fulfilled in our 
experiment, as well as in those performed by other au- 
thors Q), since for non adiabatic processes (for example, 
in fused salt) the coexistence of local solid nuclei in the 
liquid phase has been observed [HI]. 

An unexpected novel transition at T3 = 1.3 ± 0.6°C 
has been as well detected with a non negligible proba- 
bility, calling for further accurate investigation in order 
to achieve a more complete understanding of the unique 
properties of water. 
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